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ELECTRICAL CHARACTERIZATION OF 


6H CRYSTALLINE SILICON CARBIDE 
STEPHEN E. LEMPNER 


ABSTRACT 


Crystalline silicon carbide (SiC) substrates and 
epilayers, undoped as well as n- and p-doped, have been 
electrically characterized by performing Hall effect and 
resistivity measurements (van der Pauw) over the temperature 
range of approximately 85 K to 650 K (200 K to 500 K for p- 
type sample) . By fitting the measured temperature dependent 
carrier concentration data to the single activation energy 
theoretical model: 1) the activation energy for the nitrogen 
donor ranged from 0.078 eV to 0.101 eV for a doping 
concentration range of 10 17 cm' 3 to 10 18 cm' 3 , 2) the 
activation energy for the aluminum acceptor was 0.252 eV for 
a doping concentration of 4.6 x 10 18 cm' 3 . By fitting the 
measured temperature dependent carrier concentration data to 
the double activation energy level theoretical model for the 
nitrogen donor: 1) the activation energy for the hexagonal 
site was 0.056 eV and 0.093 eV corresponding to doping 
concentrations of 3.33 x 10 17 cm' 3 and 1.6 x 10 18 cm' 3 , 2) the 
activation energy for the cubic site was 0.113 and 0.126 eV 
corresponding to doping concentrations of 4.2 x 10 17 cm' 3 and 
5.4 x 10 18 cm' 3 . 
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CHAPTER 1 


INTRODUCTION 

The interest in silicon carbide (SiC) steins from 
properties that make it suitable for applications difficult 
or impossible for other semiconductors to duplicate. It is 
an adamantine refractory material with a wide bandgap. 
Further, it is chemically inert and radiation tolerant [1], 
making it suitable for use in extreme environments. Other 
outstanding properties of SiC include its high breakdown 
field E B , a high saturated drift velocity v L (v L is the 
velocity at which the electron has enough energy to emit an 
optical phonon [2], v L = V[ (8E o ) / (37rm*) ] [3], where E o is the 
optical-phonon energy and m* is the effective mass of the 
electron) , and a high thermal conductivity k (at room 
temperature and below, silicon carbide has a higher thermal 
conductivity than that of pure copper; at room temperature 
the thermal conductivity of 6H-SiC is 4.9 W/cm K as compared 
to 4 W/cm K for copper [4]). 

SiC crystallizes in several forms. This research is 
restricted to the form known as 6H-SiC. Table 8 in Appendix 
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A gives some of the parameter values for 6H-SiC at room 
temperature, along with those of Si, Ge and GaAs for 
comparison 1 . 

Several figures of merit have been developed to compare 
how the physical and electrical properties of semiconductors 
would affect their device performance. Two widely used 
figures of merit [5,6] combine the saturated drift velocity 
with the breakdown field or thermal conductivity to 
determine frequency and power limitations of the material. 

The Johnson (1965) figure of merit, Z., uses the 
breakdown field and the saturated drift velocity to arrive 
at a measure of the high frequency and high power capability 
of devices. It is written as 


Zj = 


Eb^Vl 

4 *7t 


(V/sec) : 


( 1 ) 


where Zj sets an upper limit on the combination of device 
parameters, namely, the maximum power output P m , the maximum 
operating frequency F m , and the maximum load impedance Z m 
[2,5]. 

The Keyes (1972) figure of merit, Z K , emphasizes the 
thermal conductivity limitations of devices. It is written 


1 Appendix B gives a brief history of silicon carbide 
and its nomenclature. 
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z y = k( C * Vl ) 2 (W/ sec-°C) (2) 

Here, c is the speed of light and e is the static dielectric 
constant of the material. Z K includes small size effects. 
Smaller devices have higher speed, limited by the saturated 
drift velocity, and higher device density increases power 
dissipation per unit area [2]. 

The figures of merit illustrate that the properties of 
SiC lend themselves to power, microwave, and elevated 
temperature applications; not, however, to small devices 
operating at room temperature [2], because of the large 
bandgap which would cause, for a given doping concentration, 
the width of the depletion region to be higher in SiC 
devices than in devices made of lower bandgap materials. 
This in turn would limit the theoretical size reduction 
possible. 2 

The purpose of this research was to electrically 
characterize bulk and epilayer 6H-SiC samples. The 
electrical characterization consisted of Hall voltage and 
resistivity measurements so as to determine the majority 
carrier concentrations and mobilities. The values of 
majority carrier concentration over a wide range of 
temperature allow determination of the dopant concentration 

2 Table 9 in Appendix A gives calculated values of 
Johnson's and Keyes' figures of merit for Si, Ge, GaAS and 
6H-SiC. 
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and its activation energy. Such a characterization is a 
necessary first step in both basic materials research and in 
the development and understanding of useful devices. 



CHAPTER 2 


THEORY 


This chapter presents the theory of how the carrier 
concentration and mobility values are derived from Hall 
voltage and resistivity measurements. Also presented is the 
theory of how the dopant concentration and its activation 
energy are obtained from the temperature dependence of the 
majority carrier concentration. This is followed by the 
theory of measurement accuracy and possible extraneous 
effects . 

2.1 Hall Effect 

E. H. Hall devised an experiment in 1879 [45,46] to 
determine the nature of the force acting on a material 
carrying a current in a magnetic field. He found that by 
setting a magnetic field H perpendicular to a strip of the 
material through which current I is flowing, a force 
develops perpendicular to both the current and magnetic 
field. This force is F = II X B, where 1 is the directed 
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length of the strip along the direction of current flow, and 
B is magnetic flux density. This force is due to the 
combined effect of the Lorentz force qv d X B acting on each 
charge carrier within the conductor. Here v d is the drift 
velocity of charge carriers and q is the magnitude and sign 
of the charge on the carriers. As seen in Figure 1, for the 
same direction of conventional current and the magnetic 




Figure 1 Charge Carriers in a Magnetic Field 


field, the direction of the Lorentz force is the same on 
both electrons and holes. This produces a transverse Hall 
potential 

v h = J* E H -dw = E h w 


( 3 ) 
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between the sides of the strip, where w is the width of the 
test strip and E H is the Hall electric field intensity. The 
sign of this potential determines the sign of the majority 
charge carriers. 

The Hall electric field E H builds to the point where it 
just cancels the force caused by the magnetic field. 

qE„ = g(v d x B) (4) 

For an n-type extrinsic sample, the number of charge 
carriers per unit volume, n, can be found from Hall effect 
measurements. Since v d and B are perpendicular to each 
other, E h = v d B and since 

v d = (cm/sec) (5) 


where J is the current density in the strip, the carrier 
concentration is given by 


J B 
Q e h 


( 6 ) 


or, in terms of measurable quantities 


n 


I B 
q t V H 


(cm' 3 ) 


(7) 


where t is the thickness of the test strip along the 
direction of the magnetic field. 
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The Hall coefficient R H is defined as the ratio of the 
Hall electric field to the product of current density and 
magnetic flux density [44]. It is written 


or, in terms of measurable quantities 

R = (cm 3 /C) ( 9 ) 

H IB 


Comparing Eq. (7) and Eq. (9) , the Hall coefficient R H = 1/nq. 
However, R H is commonly expressed as 


where the meaning of the Hall factor r H will become clear 
shortly . 

The signs of the Hall coefficient and the Hall voltage 
are the same; if conduction is due to electrons they are 
negative, if it is due to holes they are positive. 

Drift mobility is defined as the ratio of the drift 
velocity v d to the electric field intensity E parallel to 
the current flow, and is written 





(cm 2 /V sec) 


( 11 ) 
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Resistivity p is defined as the ratio of the electric 
field intensity, parallel to the current flowing in the 
material, to the current density. 


P 


E 

J 


( 12 ) 


Substituting for the current density J from Eg. (5) we get 


P 


E 

n q v d 


( 13 ) 


Using Eq. (10) and Eq. (11) gives 


P = 


r H P 


(n cm) ( 14 ) 


for the resistivity of the material. 

The Hall mobility n H is defined as 


^ = 


1**1 

P 


(cm 2 /V sec) 


( 15 ) 


Combining Eqs. (14) and (15) we get the relationship 
between Hall mobility and drift mobility, /x H = r H n, where 
r H is the Hall factor which is usually close to one and is 
a function of the energy-dependent relaxation time which is 
dependent on the scattering mechanism(s) involved [49]. 
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2 . 2 van der Pauw Method 

In 1958, L. J. van der Pauw published a theorem [42] 
for measuring the specific resistivity and Hall coefficient 
of flat lamella of arbitrary shape. This can be done 
without knowing the current pattern in the sample if the 
following conditions are met. 

1. The contacts are at the periphery of the 
sample. 

2. The contacts are of negligibly small area. 

3. The sample has uniform thickness and has 
no holes. 

The four needed contacts can be placed in an arbitrary 
way along the periphery of the sample, as shown in Fig. 2. 
The resistivity is calculated from 
two resistance measurements R 1 and R 2 
as defined below. When a current I 1 
is passed into A and out of B, 
measuring V 1 the potential difference 
between V 0 and V c (V, = V D - V c ) , gives 
resistance R 1 = V 1 /I 1 . In a similar 
way resistance R 2 = V 2 /I 2 is found by 
passing a current I 2 into B and out of 
C and measuring the potential difference V 2 between V A and 
v„ (V, = V, - V„) . 



Flat Sample 
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He derived the following formula to determine the 
resistivity, using these resistance measurements. 


_ t (R 1 + R 2 ) Ri ^ 

ln2 2 R 2 


( 16 ) 


where t is the thickness of the sample and the value of the 
function f(R.,/R 2 ) is obtained from the following: 

CTT * -H2 arc c° s *<| exp ( J f £)) < 17 > 

where Q = I^/Rg or its reciprocal so that Q>1, and 
f = f (VR 2 ) . 

2 . 3 Measurements 

The resistance measurements were made in two 
configurations, as shown in Figure 3, with 1=2 00 /liA. This 
covers all four sides and contacts. The measurements were 
repeated after reversing the current. 

The resistivity was calculated for each configuration 
and if the resistivities did not agree within 10%, the data 
at that temperature was discarded. When inconsistencies 
occurred, they did so only at low temperatures, possibly due 
to the probe tip not making good contact. 

The Hall voltage is determined by forcing a current 
into A and out of C and measuring the potential developed 
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voltage across B and D from when the magnetic field is 
absent to when it is present, is the Hall voltage V H . 
Equation (9) , used to calculate the Hall coefficient is 
repeated here. 

R h = -£jV h (cm 3 /C) (18) 


The current I was 1 mA and the magnetic flux density 
B was 5000 gauss. Readings were repeated after reversing 
the current. No heating effects were observed. The Hall 
magnet and power supply did not allow the magnetic field to 
be reversed. 
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2.4 Activation Energy from Carrier Concentration 

The derivation of the equations used to determine the 
activation energy from the carrier concentration data is 
presented for both cases, namely, where the dopant impurity 
is assumed to give rise to 1) a single localized energy 
level (single level model) , 2) two distinct localized energy 
levels (double level model) in the forbidden band. The 
derivation assumes non-degenerate compensated material. If 
the donor impurity site has a bound electron it is assumed 
electrically neutral regardless of the state of the 
electron, namely, the ground state or an excited state [54]. 

The following derivations will be done for n-type 
materials since the form of the equations for p-type 
material is the same, except that 8, the degeneracy factor 3 , 
is usually larger than one [54]. 

2.4.1 Single Level Model 



6 — 0.5 for donors and 8 = 2. for acceptors was used. 


3 
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where n is the electron concentration, N a * is the density of 
ionized acceptors, and N d is the total density of donors. 
The constant R is the spin degeneracy factor and depends on 
the number of ways an electron can occupy a given level. 
The reduced activation energy e d = E d /kT, and the reduced 
Fermi energy T) = E F -E c /kT. Here E d is the activation energy 
which is the same as the donor ionization energy E c -E donor , 
and F f , E donor and E c are Fermi energy level, donor energy 
level and the conduction band edge, respectively. 

Assuming 100% ionization for the minority dopant 
acceptors, equation (19) can be rewritten as 

(n+N a ) (l + P _1 exp(e d + ti) ) = N d (20) 

Since, for a non-degenerate case exp (77) = n / N c / where N c is 
the effective density of states 4 in the conduction band, 

(n + N a ) (i + -JL- exp (e d ) ) = N d (21) 

Multiplying out the left side of this equation and combining 
terms gives 


4 N c depends on the density of states effective mass of 
electrons, m^*. Appendix E gives a detailed discussion of 
the effective mass of electrons and holes in 6H-SiC. In the 
data reduction m^* = 0.915 m o and the density of states 
effective mass of holes iru* = 1.0 m . 

an o 
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n 2 [-j—exp(e d ) ] + 

Dividing through by n 
quadratic formula gives 

- [l + -j^-exp(e d ) ] -/T[l + 

1 _ P Nc 

n 


n [1 + - 
2 and 


^-exp(e d )] + [N a -N d ] 
solving for 1/n 


= 0 ( 22 ) 
using the 


Therefore, 


n=- 


2 (N d -N a ) 


N. 


N, 


[l + - i r£-exp(e d ) ] +v /T[i + _£-exp(e d ) ] 2 + 


pW, 




4 (iV d -rt a ) 
P^c 


exp (e d ) ) 


(24) 


which is the equation used in the least squares analysis 
(Appendix C) . E d , N d and N a are the parameters solved for by 
an iterative method. 

2.4.2 Double Level Model 


The general form of the equation for compensated 
material is given by 


*a = 


t-md 

Z^j-i 


Ns 


1 + Pj 1 exp(e d ^ + ti) 


n + 


(25) 
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which gives the following cubic equation in n for the double 
level analysis (m=2) . 


exp(e d +e d ) exp(e d ) exp(e d ) exp (e d +e d ) 

n 3 [ . ] +n 2 [ — „ ■ — t „ - + n~ n — ] + 


PiP 2 ^; 

N-N, 


Pi 

N-N, 


P 2 ^c 


PlP 2 ^C 


(26) 


n[1+ p i ^ 2 exp(e d 1 ) + -^- iex P< e d 2 )i + t^ a --^dr^d 2 ] =° 


The Italian mathematician Girolamo Cardano [55] found 
the closed form solution for a cubic of the form 

x 3 +px+ g=0 (27) 


to be 


x = A+B 

— — A+I3 . A-B /_2 (28) 

" 2 2 * 


where 


A 

£ 


' 2 + \ 

g 2 , P 3 
4 27 

"2~\ 

g 2 + P 3 
4 27 


(29) 


Equation (26) can be expressed as 


n- 


a b 

— L + — + c = 0 
n 2 n 


(30) 


and it can be made to resemble equation (27) by substituting 
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_1 

n 



into the above cubic equation. 
This gives 

a 2 


(31) 


(32) 


and 


Q = 


2 a 3 
27 


ab 

3 


+ c 


(33) 


Substituting equations (32) and (33) in equation (29) and 
the resulting values of A and B in equation (28) , we can 
solve for n from equation (31) . 

The values: 1) N d1 , N d2 , the doping concentrations 
corresponding to the two ionization (activation) energy 
levels E d1 and E d2 respectively, 2) E d1 , E^, and 3) N a are 
obtained iteratively from a least-squares fit of the 
temperature-dependent carrier concentration data. 

2.5 Activation Energy from n(T) and p(T) at low 
Temperatures . 


Equation (21) for n is repeated here 
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(n + N a ) (l + ^ r exp(e d ) ) = N d 


( 34 ) 


With N a = 0 and at low temperatures, where kT << E d , equation 
(34) can be approximated as [56] 


-E a 


n = yjWFd exp( ^| ) 


( 35 ) 


Hence, E d can be found from the slope of ln(n) vs 1/T 
at low temperatures . 

Equation (13) for the resistivity is repeated here 


p = £ m __g_, 
J n g v d 


( 36 ) 


Substituting Eq. (35) into Eq. (36) gives [46] 


P = 


B eXP( 2^ } 

/POT 


( 37 ) 


Again, E d can be found from the slope of ln(p) vs 1/T at low 
temperatures . 

In equation (35) and (37) , the temperature dependence 
of the exponential factor completely masks the temperature 
dependence of the pre-exponential factor. Hence, the pre- 
exponential factor can be treated as a constant with respect 
to the temperature. 
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2 . 6 Measurement Accuracy 

2.6.1 Error Due to Contact Size 

A systematic error is introduced in the measurements 
when the sizes of the contacts are not negligible in 
comparison with the dimensions of the sample. By estimating 
the error [47], a correction factor can be applied to the 
results. For the samples used in this research these 
corrections were not a factor. 

2.6.2 Calculation of Error in the Final Result 

During a measurement, the voltmeter (Fluke 8520A) had 
the capability of being programmed to take multiple readings 
and do a statistical analysis. The t-distribution was used 
to obtain the uncertainty in the sample standard deviation 
s . The confidence limits on the measurements x are ±t ,s„, 
where a represents the confidence level (i.e. 90%, 95%, 99%) 
and <p = n-1 is the degree of freedom in the measurement, 
where n is the number of measurements. There were 61 
measurements so for the 95% confidence level a = .025 and 

025, 60 = 2 ' 

The error propagation was calculated from the initial 
measurements to the final calculated result. Since the 
variance on the primary inputs, the thickness and the 
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voltage, was known the variance on the results could be 
calculated, using the following relationship from [48]: 


V „ 


(l£) 2 r + ( i £ )2 ^ 

at c v v 


dv 


( 38 ) 


As an example, to find the confidence limits on the 
carrier concentration n, where n is defined by 


n = 


( q r h 


( 39 ) 


and a stands for r H /q. Then 


dn 


dR 


H 


R H 


n 

Ru 


( 40 ) 


so that from equation (38) 


v„ - (-^-) 2 V. 


R 


( 41 ) 


To calculate V RH , the variance in R,,, we first express R H as 


R u = 


B I 


= a v u t 


( 42 ) 


where a stands for (1/BI) . 
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dv H 


= a t 


£h 

V H 


( 43 ) 


dR H 

dt 


a v H 


t 


( 44 ) 


Then, using Eq. (38) gives us the following expression for 
V RH , the variance in R H : 


= A 2 


V„ 


(J T )2 


( 45 ) 


The final result for the variance on the carrier 
concentration is obtained when equation (45) is substituted 
into equation (41) . 




( 46 ) 


Then, the standard deviation in n is given by 


s 


n 



( 47 ) 


and the limits on n for 95% confidence level are n±2s . 

n 

The variance of the voltage, V yH was calculated from 
voltmeter readings; V t was determined as follows. The 
thickness of the 4/um epilayers was accurate within ±l/xm. 
Although this was arrived at through many (hundreds of) 
observations, it was justified in letting <p equal 60. With 
this decision, ±t a<# s t = ±2s t = ± 1 ^m so that s t =0.5 and 
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V =0.25. This was used in the calculation of the 
confidence limits on the resistivity p, carrier 
concentration n, and Hall mobility /x H - Figures 4 through 7 
show n plotted against reciprocal temperature. At each 
temperature, the variance in n is indicated by an error bar. 

The voltmeter and electrometer were considered 
accurate. The voltmeter error was approximately ±0.012% on 
the lOOmV scale, and the electrometer error was ±0.5% on the 
10' 3 A scale. 

2.7 Extraneous Effects 

The measurements can be influenced by thermoelectric 
and thermomagnetic effects besides the galvanomagnetic 
effect we are trying to measure. 

The thermoelectric effects are the Seebeck and Peltier 
effects. The Seebeck effect is the effect produced by a 
hot-point probe where the mobile charge carriers move away 
from the higher temperature area. Both the type and 
quantity of charge can be measured as in the Hall effect. 
The voltage produced by the Seebeck effect is V s = -a Q AT, 
where a 0 is the thermoelectric coefficient and AT is the 
temperature difference. This voltage can be measured and 
accounted for while making the Hall voltage and resistivity 
voltage measurements. 
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Figure 4 Carrier concentration as a function of temperature 
for the undoped Cree substrate B04365-5B with error bars. 


Log Carrier Cone vs Reciprocal Temp 

Undoped Sample B0428-2H 



1 ** 4 ®«7*e io ii is is 

Reciprocal Temperature (lOOO/T) K ** 1 

Figure 5 Carrier concentration as a function of temperature 
for the undoped Cree substrate B0428-2H with error bars. 
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Figure 6 Carrier concentration as a function of temperature 
for the undoped NASA grown epilayer 1312-3 with error bars. 



Figure 7 Carrier concentration as a function of temperature 
for the aluminum doped NASA grown epilayer 1270-2 with error 
bars. 
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The Peltier effect results from the different amounts 
of heat that are carried by the current in different 
materials. If two materials with different Peltier co- 
efficients are in intimate contact, the junction will either 
get warmer or cooler depending on the direction of the 
current and the difference between their Peltier 
coefficients. The magnitude of the current and the 
magnitude of the difference between the Peltier coefficients 
determines the magnitude of the heat current. 

The net rate of increase of heat energy at the 
junction, per unit area of the junction, is given by 

= J 0i ” J Qo = ( *1 ~ n 2> J (48) 

Here the J Q 's are the heat currents, the w's are the Peltier 
coefficients of the two materials (it, is for the material 
the electric current is entering and n 2 is for the material 
that the current is leaving) , and J is the electric current 
density. 

If any temperature gradient in the sample is caused by 
this Peltier effect, then the resulting Seebeck voltage due 
to such a temperature gradient can be eliminated by 
reversing the sample current and taking the average of the 
resistive voltage readings across the sample. 

The magnetoresistance is eliminated by taking a reading 
and reversing the direction of B and taking another reading, 
then averaging the readings. The Nernst and Rigi-Leduc 
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effects can be eliminated in a similar way by reversing both 
B and J and averaging the readings. However, the 
Ettingshausen effect cannot be eliminated by reversing and 
averaging the measured voltages. 

In short, all of the extraneous effects (with the 
exception of the Ettingshausen effect [44]) can be 
eliminated by reversing current and/or magnetic field as 
mentioned above. 

Since in this research the magnetic field could not be 
reversed, the size of the effects need to be determined to 
see the magnitude of their influence on the measurements. 

The references [3,9,46 and 49] were heavily relied upon 
for the equations used to calculate the magnitude of the 
influence of the thermoelectric and thermomagnetic effects. 

In Table 1, the calculated values of the various 
extraneous voltages and temperature differences developed as 
a result of the thermoelectric and thermomagnetic effects 
for the 6H Sic undoped epilayer sample 1312-3 are presented 
at three different temperatures. As seen from this table, 
the magnitudes of the induced voltages and temperature 
differences are negligible. 
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Table 1 Calculated Thermo- (electric and magnetic) Effects 



measurement 
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CHAPTER 3 


SAMPLE PREPARATION 


This chapter describes: 1) the 6H-SiC samples used in 
this research, 2) contact metalization, and 3) the design of 
a reactive ion etch mask to define the sample geometry. 

3.1 Growth of 6H-Sic Epilayers 

The 6H-SiC n- and p-type epilayers were grown by 
chemical vapor deposition (CVD) , in the Engine Sensor 
Technology Branch of the NASA Lewis Research Center 5 . The 
4/xm epilayers 6 were grown on 6H-SiC substrates, 
approximately 5.5 mm X 5.5 mm which were 3° off the 
crystallographic (0001) plane. The substrates were grown by 
the modified Lely method [40] by Cree Inc. 


5 National Aeronautics and Space Administration, Lewis 
Research Center, 2100 Brookpark Rd., Brookpark OH 44135. 

6 The chemical vapor deposition (CVD) grown epilayers 
were made available to me for this research courtesy of Mr. 
J. Anthony Powell. 
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The test epilayers were always isolated from the 
substrate by a 4 jum thick buffer layer as shown in Figure 8. 
In order to obtain planar test epilayers, sample edges were 
trimmed with a watering saw, because the buffer and test 
epilayers grow all around the substrate during the growth 
process. Reactive Ion Etch (RIE) patterning, if available, 
is a more convenient and more accurate method of defining 
sample geometry. 



3 . 2 Contacts 

A series of contacting schemes was tried using the 
materials Al, Au,Mo,Ni, Si,Ta, and Ti, either singly or in 
combinations. These materials were deposited either by 
evaporation or sputtering. Several of these schemes looked 
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promising. What was settled on for this research was 
aluminum for the p-type samples and nickel for the undoped 
and n-type samples. The contacts were Rapid Thermal 
Annealed (RTA) in forming gas (4% Hydrogen in Nitrogen) 
until they were determined to be ohmic. 7 Generally, the 
anneal temperature was 925°C and the time was varied, from 
a few seconds to 1 or 2 minutes, to accomplish the anneal. 
It was observed that the n-type samples became a deeper 
green [1,8,41] after the RTA process. 

It was noticed (during the trial and error process of 
finding the proper time/ temperature combination for the RTA 
contacts) that if, for example, the aluminum contacts were 
overheated and the sample was etched to remove all remnants 
of the metalization, and if the new metalization overlapped 
areas where the old contacts had been, the new contacts were 
either already ohmic or required a very mild anneal to 
become ohmic. Even with these "mild", for Sic, annealing 
conditions, contrary to what is understood about the 
difficulty in doping Sic, it appears that doping had 
occurred, even if extending to only a few atomic layers. 


7 I would like to express my thanks to Carl Salupo for 
his assistance in the metalization and annealing processes. 
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3 . 3 Mask 

Although not used in this research, a mask set was 
designed for Reactive Ion Etching (RIE) of van der Pauw 
patterns and Hall bars on the Sic samples for electrical 
characterization studies. The mask had two Hall bars and two 
van der Pauw patterns within a 2mm X 2mm square as shown in 


TO ip Jl wn TO ip Jl fl TO ip 
ip TO fl aTl ip TO Jl Jl ip TO 


Figure 9 Mesa Mask Layout 
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The mask set was designed to reduce the measurement 
errors caused by contact size and placement (i.e. the 
contacts not being infinitesimally small and not being on 
the periphery of the sample [42]). 

Two primary sources used for the design criteria of the 
Hall bars and van der Pauw patterns were Laboratory Notes on 
Electrical and Galvanomaanetic Measurements [43] for the 
orthogonal pin-wheel van der Pauw pattern dimensioning, and 
The American Society for Testing and Materials (ASTM) 
Designation: F76-86 [44] for the six contact Hall bar 

dimensioning. 

There are two objectives that need to be balanced in 
the design, namely, the dimension ratios to reduce the 
errors caused by contact placement and the dimension ratios 
to account for the range of material resistivities that will 
be encountered. For example 
as seen in Figure 10, the 
pin-wheel leg length (A) 
needs to be long compared 
to the sample area 
dimension (S) , and the leg 
width (B) needs to be small 
compared to this same area 
dimension (S) . As either 
ratio moves in the 
indicated direction the 


, in the van der Pauw pattern. 



Figure 10 

Pattern 


RIE Pin-wheel 
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measurement error decreases [43], but as the resistivity of 
the material increases these ratios need to move in the 
opposite direction to reduce heating effects. 

The actual dimensions of the test patterns (Figures 11, 
12) took into consideration the above-mentioned dimension 
ratios as well as the minimum area needed for the 
wirebonding pad size and the capabilities of the clean room 
photolithography environment which determines the minimum 
attainable line width. 

We next consider how the upper limit of sustained 
current flow through a sample relates to pattern dimensions. 
What should the input power limitation be when measuring Sic 
epilayers so that the temperature rise in the specimen is 
limited to 0.5°C? Note that ASTM allows a ±1°C. 

We know the power (P) dissipated by the sample is 



where I is the current through the sample, p is the 
resistivity of the specimen, L 1 is the length of the current 
carrying leg, and A 1 is the cross sectional area of current 
flow. We also know the heat flow through a sample per unit 
time is given by 


P = 


x A 2 


AT 

L 2 


( 50 ) 



Figure 11 Dimensioned Hall Bar 


where k is the thermal conductivity, Aj is the area for heat 
conduction in the sample, AT is the temperature difference, 
and L 2 is the length of the heat flow path, in this case 
t/2 . 

Combining Eq. (49) and Eq. (50) gives 


J 2 p 


K A 1 A 2 A T 
2 


(51) 







Figure 12 Dimensioned van der Pauw Pattern 

for the van der Pauw configuration. The factor 2 in the 
denominator comes from the two current carrying legs of the 
sample. The leg dimensions will be the limiting factor in 
the design. 

The areas are A 1 = (w t) and A 2 = 2 (w L t ) , where the 
factor 2 on the right hand side of A 2 is because the heat 
will be dissipated through both the top and bottom of the 
legs. Here w is the width of the current carrying leg (35 
microns for the van der Pauw and 70 microns for the Hall bar 
patterns, as shown in figures 11 and 12) and t is the 
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epilayer thickness. Substituting for A,, A 2 and L 2 in 
equation (51) gives, for the van der Pauw pattern 

J 2 p = 2k w 2 AT (52) 

and for the Hall bar gives 

J 2 p = 4k w 2 AT (53) 

Table 2 below gives values of I and p satisfying these 
equations for AT = 0.5 K. 


Table 2 Current and Resistivity Values for AT = 0.5 K from 
Equations (52) and (53) 


CURRENT 

RESISTIVITY UPPER LIMIT (n-cm) 

ma 

van der Pauw 
w = 35 nn 

70 K 300 K 1000 K 

Hall 

w = 70 ixm 

70 K 300 K 1000 K 

200 

13E3 

1500 

330 

10E4 

12E3 

2600 

500 

2100 

240 

53 

16E3 

1900 

430 

1000 

520 

60 

13 

4200 

480 

100 
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Table 3 below gives, at each temperature, of the five 
samples used in this research, the sample which had the 
highest resistivity and value of the resistivity. 

Table 3 Highest Resistivity Encountered 


Temperature 

p 

Sample 

K 

n-cm 

Number 

82 

132.7 

B0428-2H 

300 

2.01 

1270-2 

650 

0.82 

1270-2 


Comparing tables 2 and 3, it is easy to see that for 
the samples used in this research, the dimensions defined by 
the RIE patterns would cause the temperature rise to be 
significantly smaller than 0.5 K for the measuring current 
I in the range of 200/iA to 1mA. 

The mask set was designed and the masks were received. 
However, the RIE was unavailable for this research. So, the 
samples that were tested with a grown epilayer were trimmed 
with a wafering saw to eliminate the effects of growth 
encapsulation mentioned earlier. The samples were 5.5mm X 
5 . 5mm square and the contacts were placed in the corners of 
these samples. For these dimensions the temperature rise AT 
caused by the measuring current in the range of 200/xA to 1mA 
should be orders of magnitude smaller than that calculated 
for the RIE patterns. 
















CHAPTER 4 


APPARATUS 


4.1 Hall Setup 

The magnet was a Varian Associates 6" Hall magnet 
(V3701F) and its power supply (V-FR2603) was equipped with 
a magnetic field regulator. 

The magnetic field was monitored with a rotating-coil 
gaussmeter 8 . This gaussmeter was calibrated to read 2.94kG 
using a standard magnet 9 of 2938 gauss. The Hall 
measurements were done at 5kG. 

For temperature-dependent Hall measurements, the 
temperature was measured with a type J thermocouple 10 


8 Rawson Electrical Instrument Co., Cambridge, Mass., 
type 820, number 17270. 

9 Ibid 

10 Type J thermocouples are iron-constantan . 
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in conjunction with a Doric Trendicator 400A. The 
thermocouple had a temperature range from -200 to +933 °C, 
and the Doric Trendicator gave a one degree Celsius 
resolution for the temperature reading. 

4.2 Cryostat and Sample Jig 

It was necessary to design and build a cryostat and 
sample platform for temperature-dependent electrical 
measurements on the test samples. For the design, it was 
deemed necessary that the magnetic metal content be 
minimized. Most of the design was accomplished with 
available off-the-shelf materials. 

The cryostat was fabricated from a firebrick 11 , which 
was cut in half and machined to 

1. provide areas for the sample platform, gas 
torch 12 , sample test leads and thermocouples; 

2. fit between the pole pieces of the Hall 
magnet, which had a 2" clearance. 

Appendix D shows detailed schematics of the cryostat 
construction. 


n Nock Fire Brick Co., series K-23. 

12 GTE, serpentine gas/air heat torch part 014923. 
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The test jig was a miniature probing station comprised 
of the sample platform and probes as shown in Figure 13. 
The sample platform was 
fabricated from Macor 13 , a 
machinable glass-ceramic. The 
probes were made from 
nonmagnetic stainless steel 
washers, with slotted center; 
spot welded to these was a short 
section of nonmagnetic stainless 
steel wire to which standard 
tungsten probing tips were spot 
welded. The probes provided three degrees of freedom. 

4.3 Temperature Measurement £ Control 

There were two thermocouples, one in contact with the 
bottom of the sample and the other approximately 3mm above 
the sample platform. The two thermocouples were allowed to 
stabilize before a measurement was taken. The samples were 
extremely sensitive to temperature variations. For 
meaningful readings it was necessary that the temperature 
remained constant within «±1°C during a measurement. 

The low temperature measurements were controlled by 
adjusting the flow of nitrogen from a reservoir, which 



Figure 13 Sample 

Platform and Probe 


13 A product of Corning Glass Works, Corning, N.Y. 
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approximated an infinite source, through a coil immersed in 
liquid nitrogen. 

For high temperature measurements, the voltage across 
the gas torch was controlled while nitrogen flow was kept 
constant. 

4.4 Electrical Measurement 

The remainder of the measurement apparatus consisted of 
a Fluke 852 OA digital multimeter, used to measure the 
developed voltage, a Keithly Instruments 225 current source 
and a Keithly 616 digital electrometer, to measure current 
through the sample. These were connected to the sample 
through a switching box used to make the van der Pauw 
connections as shown in Figure 14. The square in the upper 



lower portion of Figure 14 van der Pauw Connections 
the figure 
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represent the possible switch positions; positions 1 and 2 
are for resistivity measurements, position 3 is for the Hall 
measurement. 



CHAPTER 5 


Results and Discussion 

In this chapter is presented: 1) a very brief 

historical background on the rise of silicon carbide to the 
point of being a viable semiconductor material for the 
fabrication of useful devices, 2) activation energies 
obtained by previous researchers using both the single and 
double activation energy level models, and 3) results of 
this research along with discussion, and 4) recommendations 
for future work. 

5 . 1 Review of Past Work 

The following historical information has been taken 
from [7], [8], and [9]. 

The semiconductor nature of silicon carbide was first 
investigated early in the twentieth century. Until 1955, 
the crystals used in investigations were results of chance 
occurrences from the commercial (Acheson) process of 
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manufacturing silicon carbide as an abrasive. Generally, 
neither the polytype 14 nor the dopants were known for sure. 

There was a thorough study made by Busch and Labhart in 
1946 on Acheson crystals, and by Lely and Kroger in 1958 on 
crystals grown and doped by the process Lely published in 
1955. 

The Lely method was an attempt to produce purer 
crystals of known doping concentration on a laboratory 
scale. This was the impetus for the flurry of silicon 
carbide investigations that occurred during the sixties and 
the early seventies. 

What finally raised silicon carbide to the place of 
being a viable material for device fabrication was the 
publication of two papers, [10] in 1981 and [1] in 1983, on 
what is called the modified Lely method. 15 This brought 
about the present renewed interest in silicon carbide. 

Results of Hall measurements on silicon carbide have 
generally been based on modeling Sic as either a single 
activation energy system or as a double activation energy 
system. These results are presented next. 


14 See Appendix B. 

15 See Appendix B. 
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5.1.1 Single Activation Energy Impurity Level System 

S. H. Hagen and C. J. Kapteyns [16] determined the 
ionization energy of nitrogen donors in 6H-SiC that ranged 
from 0.081 to 0.095 eV depending on the dopant 
concentration . 

Gerhard Pensl et al. [17] found an activation energy of 
0.070 eV for heavily doped (low 10 +19 cm' 3 ) nitrogen implanted 
6H-SiC samples. 

T. Tachibana et al [18] found a donor ionization energy 
of 0.084 eV for undoped 6H-SiC. 

Lomkina et al [15] found an activation energy of 0.095 
eV for nitrogen donors in 6H-SiC with a dopant concentration 
of -1.3X10* 17 cm' 3 . For aluminum acceptors, they found an 
activation energy of 0.24 eV for an impurity concentration 
of -5.0X10* 17 cm' 3 . 

Most authors reported values for the ionization energy 
of nitrogen ranging from 0.070 to 0.095 eV 
[7,9,15,18,19,20,21,22,23,24,25] and for aluminum from 0.190 
to 0.280 eV [7,8,9,12,15,23,24,25,26,27,28]. 

5.1.2 Double Activation Energy Impurity Level System 

The existence of more than one activation energy (Fig. 
15) , for dopants in 6H-SiC was explained in [11] and [29] by 
the fact that there are three inequivalent lattice sites for 
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impurity atoms. The number of inequivalent lattice sites is 
polytype dependent [11,23]; the larger the unit cell the 
more inequivalent sites present. The inequivalence 



manifests itself beyond the nearest neighbors [24] in the 
lattice stacking sequence. 

The three inequivalent sites in 6H-SiC are categorized 
as a hexagonal site and two cubic sites. There is a larger 
binding energy difference between the hexagonal site and the 
two cubic sites (approximately 50-60 meV) 
[21,24,30,31,32,33] than between the two cubic sites 
themselves (approximately 5 meV) [21,30]. This justifies 
the use of the double activation energy model in spite of 
the fact that there exist three distinct activation energy 
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levels. Hall measurements have only resolved between the 
hexagonal and the two cubic sites, but has not resolved the 
levels between the cubic sites themselves. 

These inequivalent energy levels do not conform, as 
many authors believe, to the unaltered hydrogen atom model 
[16,19,24,32]. A reason for this is that the hydrogen model 
does not distinguish between the inequivalent lattice sites 
[24], There is also the difference caused by whether an 
impurity atom replaces silicon, as aluminum [29] (E a =0.27eV) 
usually does, or replaces carbon, as nitrogen [19,29] 
(E d =0.15 eV) and boron [19] (E a =0.39 eV) do. 

Table 4 shows, for nitrogen and aluminum in 6H-SiC, the 
activation energies that are obtained from electrical 
measurements using the two activation energy level model, 
along with the activation energies obtained optically. 

5.2 Results and Discussion 

There were five samples evaluated, two were virgin Cree 
substrates, and three were NASA-grown epilayers, one 
undoped, one n-type, and one p-type. 

The n-type NASA epilayer sample 1314-3 was grown with 
2 atoms of nitrogen per million of the host (Si and C) atoms 
in the ambient atmosphere of the growth chamber. 
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Table 4 Two Level Activation Energies 


Nitrogen (eV) 


Aluminum (eV) 


0.075 


0.100 


0.135 


0.155 



0.101 (0.158,0.163) 


0.080 


0.096 


0.100 


0.130 


0.142 


0.150 


0.081 (0.138,0.142) 


0.080 


0.086 


0.063 


0.094 


0.085 


0.120 


0.125 


0.120 


0.118 


0.119 


Method Ref 


optical 24 



0.249 optical 23 


optical 


electrical 


electrical 32 


^dh 1S hhe hexagonal site donor activation energy 
E dc is the cubic site donor activation energy. 

E ah the hexagonal site acceptor activation energy. 
E ac is the cubic site acceptor activation energy. 


electrical 33 


optical 30 


electrical 30 


electrical 


electrical 


The p-type NASA epilayer sample 1270-2 was grown with 
2 aluminum atoms per million of the host (Si and C) atoms in 
the ambient atmosphere of the growth chamber . 
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All three undoped (at least, not intentionally doped) 
samples were found to be n-type from Hall measurements. The 
n-doped and p-doped epilayers were found to be n-type and p- 
type respectively. 

Table 5 gives room temperature values of the electrical 
resistivity p, the majority carrier Hall mobility m h and the 
majority carrier concentration for the test samples. 

Table 5 Room Temperature p, and carrier concentration 
for the test samples. 



5 1270-2 


SAMPLE 

Descrip. 

SAMPLE 

TYPE 

P 

n-cm 

cm 2 /V-s 

Carrier 

Cone. 

cm' 3 

Undoped 

Substrate 

n-TYPE 

.077 

161.4 

5 . 02el7 

Undoped 

Substrate 

n-TYPE 

.088 

157.6 

4 . 50el7 

Undoped 

Epilayer 

n-TYPE 

.250 

235.8 

1. 06el7 

N 2 -2ppm 

Epilayer 

n-TYPE 

.042 

275.9 

5 . 41el7 

Al-2ppm 

Epilayer 

p-TYPE 

.250 

57.3 

5 . 40el6 
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5.2.1 Single Activation Energy Level Analysis 

Figures 16, 18, 2 0 and 22 show the carrier 
concentration versus T' 1 for the two undoped Cree substrates, 
and the NASA grown undoped and NASA grown p-doped epilayers 
respectively. The solid line represents curve-fitting using 
equation (24), valid for the single level model. The fit is 
reasonably good for the two Cree undoped substrates B0436-5B 
(Figure 16) and B0428-2H (Figure 18). For the NASA grown 
undoped epilayer 1312-3 (Figure 20) the fit is reasonably 
good for low temperatures but the data diverge from the 
theory above room temperature. For the NASA grown p-doped 
epilayer 1270-2 (Figure 22) the fit is not very good above 
500 K and below 200 K. 

There was not enough data taken on the NASA grown n- 
type epilayer sample 1314-3 to do a meaningful analysis. 

Table 6 gives the activation energy and the dopant 
concentrations N. and N. obtained from curve-fitting the 
temperature— dependent carrier concentration data to equation 
(24) , valid for the single level model. For the activation 
energy, E d of -nitrogen, this research found a range of 0.078 
eV to 0.101 eV as compared to the range of 0.07 eV to 0.095 
eV reported in the literature. For aluminum, an activation 
energy E a of 0.252 eV was found which is well within the 
range of 0.190 eV to 0.280 eV reported in the literature. 
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Table 6 Calculated Activation Energy and Impurity 
Concentrations from the Single Level Model. 


Activation Energy 
from low temp 



log n (T) 


0.202 eV 0.247 eV 


0.196 eV 0.225 eV 


0.097 eV 0.113 eV 


0.267 eV 0.243 eV 
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Log Carrier Cone vs Reciprocal Temp 

Undoped Sample B0436-6B 



Reciprocal Temperature ( 1 0OO/T) K" 1 


Figure 16 Carrier concentration as a function of 
temperature for the undoped Cree substrate B0436-5B. 



Log Resistivity vs Reciprocal Temp 

Undoped Sample B0436-6B 



* i 1 1 1 1 1 f 1" 1" " " '1 " 


1 o * 

r A •! 

i a : 

I a ; 

A 

1 

lO 1 

r ^ n 

\ A | 

' A ‘ 

I 

1 

Ax 

4 

| 

G 1 O * 

Q. 

A 

1 1 
‘ . A 1 

1 

io- 1 

A - A 

r A a A ^ a A } 

1 1 1 1 1 I 1 1 1 1 1 

1 

1 

2 a a a a 7 a a to 11 12 ia 

Reciprocal Temperature (IOOO/T) K _1 


Figure 17 Resistivity as a function of temperature for the 
undoped Cree substrate B0436-5B. 
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Log Carrier Cone vs Reciprocal Temp 

Undoped Sample BO A 2 8 — 2 H 



Reciprocal Temperature ( 1 0OO/T) K" 1 


Figure 18 Carrier concentration as a function of 
temperature for the undoped Cree substrate B0428-2H. 


Log Resistivity vs Reciprocal Temp 

Undoped Sample BQ428-2H 



Reciprocal Temperature { 1 0OO/T) K -1 


Figure 19 Resistivity as a function of temperature for the 
undoped Cree substrate B0428-2H. 
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Figure 20 Carrier concentration as a function of 
temperature for the undoped NASA grown epilayer 1312-3. 



Figure 21 Resistivity as a function of temperature for the 
undoped NASA grown epilayer 1312-3. 





Log Carrier Cone vs Reciprocal Temp 

2 ppm Al In Ambient Sample 1270 — 2 
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Figure 22 Carrier concentration as a function of 
temperature for the aluminum doped NASA qrown 
epilayer 1270-2. 
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Figure 23 Resistivity as a function of temperature for the 
aluminum doped NASA grown epilayer 1270-2. 
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5.2.2 Results of Activation Energy From Low 
Temperature n(T) and p(T) 

The plots of log p versus 1/T for the two Cree 
substrate samples B0436-5B (Figure 17) , B0428-2H (Figure 19) 
and the NASA grown undoped epilayer 1312-3 (Figure 21) look 
normal. The plot of log p versus T' 1 for the NASA grown p- 
type epilayer 1270-2 does not appear as well behaved as the 
other three samples at low temperatures. 

Table 6 also gives values of activation energies found 
from the slope of 1) low temperature log n(T) versus l/T 
(Eq. (35)) and 2) low temperature log p(T) versus 1/T (Eq. 
(37)). As expected from equations (35) and (37), the 
activation energies obtained from low temperature n(T) and 
p (T) are in good agreement with each other. However, only 
for the p-type epilayer does the activation energy found 
this way agree with that found from fitting Eq. (24) to the 
carrier concentration data over the entire temperature 
range. 

5.2.3 Double Activation Energy Level Analysis 

In chapter 2, Eq. (26) relating the temperature 
dependent carrier concentration to the parameters of the 
double activation energy level was derived without 
accounting for excited states and valley orbital splitting. 
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As a result, fitting this equation to the measured 
temperature— dependence of carrier concentration did not 
yield meaningful results. 

Figures 24, 25, and 26 show the carrier concentration 
versus T' 1 for the two undoped Cree substrates (B0436-5B and 
B0428-2H) , and the NASA-grown undoped epilayer (1312-3) 
respectively. The solid line represents curve-fitting 
using a variation of equation (26) that takes into account 
excited states and valley orbital splitting [57]. 

Table 7 shows the results of this curve-fitting. The 
activation energies E h and E c for the hexagonal and cubic 
sites found in this table are generally comparable to those 
in Table 4 obtained from the literature, except for those 
for the sample B0428-2H. The published range of E h values 
is from 0.063 to 0.1 eV and E c values is from 0.118 to 0.155 
eV. As seen in Table 7, the values of E h and E c for the n- 
type sample B0428-2H are considerably lower than the 
published values, whereas for the n-type samples B0436-5B 
and 1312-3 they are comparable to the published values. 
This discrepancy is, to an extent, understandable if it is 
taken into consideration how difficult it is to do five- 
parameter curve fitting. 

Unlike the results reported by W. Suttrop et al [30], 
this research did not find N h :N c - 1:2, where N stands for 
donors or acceptors and the subscripts h and c stand for 
hexagonal and cubic sites. In the published literature, the 
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Log Carrier Cone vs Reciprocal Temp 

Undoped Sample B0436-5B 



Figure 24 Carrier Concentration as a function of 
temperature (double level model, sample B0436-5B) . 


Log Carrier Cone vs Reciprocal Temp 

Undoped Sample B0428-2H 



25 Carrier Concentration as a function of 
temperature (double level model, sample B0428-2H) . 
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Figure 26 Carrier Concentration as a function of 
temperature (double level model, sample 1312—3) . 


Table 7 Calculated Activation Energy and Impurity 


Concentrations from the Double Level Model. 


























60 


ratio N c /N h is found to vary between 1.4 and 77, with the 
majority of values of this ratio within 210.5. This 
research found this ratio to be 3.4 for the sample B0436-5B 

and 1.25 for the sample 1312-3. 

From the published values of E h and E c for nitrogen in 
6H-SiC, the difference E c -E h ranges from 30 meV to 57 meV. 
This research found this energy difference to be 33 and 57 
meV for the samples B0436-5B and 1312-3 respectively. 

It was not possible to fit the double level model to 
the temperature-dependent carrier concentration data of the 
p-type sample 1270-2 due to the limited reliable data. Most 
researchers do not find a double level result using 
electrical characterization for p-type material. 

5.2.4 Results of M h ( t ) 

The mobility data is presented in Figure 27. The top 
three curves are for undoped (n— type) samples and the bottom 
curve is for the p-type sample. All four samples show a 
decrease in mobility as the temperature is increased. For 
the undoped samples, the room temperature mobility (also 
given in Table 5) of electrons decreases with increasing 
total (N d + N a ) impurity concentration found from single 
level analysis. The room temperature hole mobility for - 
10 18 /cm 3 p— doped sample is about 1/3 of the room temperature 
electron mobility for comparably n-doped samples. 
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Log Mobility vs Reciprocal Temp 

+ BO4S0-OB A B0428-2H O 1J12-9 O 1270-2 



Figure 27 Mobility as a function of temperature for the 
four test samples. 


5 . 3 Summary 

Both the single and double activation energy level 
models have been used in this research to characterize 6H- 
SiC substrates and epilayers by fitting the models with the 
temperature dependent carrier concentration. The activation 
energies found from both these models are, in general, in 
fair agreement with those found in the literature. 
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5.4 Recommendations for Future Work 

A. More research is needed to firmly establish the 
double activation energy model in terms of 1) 
finding the accurate values of the cubic and 
hexagonal activation energy levels, 2) finding the 
dopant concentration for each activation energy 
level and finding the dopant dependence of the 
activation energies, for both p- and n-type SiC. 

B. The accuracy and efficiency of the measurements 
can be further improved by doing the following: 

1. ) using an automated Hall rig for data 

acquisition, 

2. ) using photolithographically mesa etched 

epilayers for accurately defining the 

sample geometry, 

3. ) using ultrasonically wire-bonded 

contacts, 

4. ) using non-magnetic thermocouples, and 

5. ) using a magnetic power supply that allows 

reversing the direction of the magnetic 


field. 
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APPENDIX A 


Tables 


Table 8 Parameter Values 1 (at room temperature) 


Parameters 

Si 

Ge 

GaAs 

6H-SiC 

Bandgap Eg (eV) 

1.12 

0.66 

1.42 

2.93 





[50] 

Thermal Conductivity 

1.5 

0.6 

0.46 

4.9 

k (W/cm K) 




[4] 

Optical-Phonon Energy 

0.063 

0.037 

0.035 

0.1205 

E 0 (eV) 




[6] 

Saturated Drift 

1X10 7 

6X10 6 

2X10 7 

2X10 7 

Velocity v L (cm/sec) 

[5] 

[2] 

[6] 

[6] 

Breakdown Field 

3X10 5 

1X10 5 

4X10 5 

2X10 7 

E b (V/cm) 




[51] 

Dielectric Constant 

11.9 

16.0 

13.1 

j.9.66 

(Static) e 




[52] 


1 All of the parameter values were taken from reference [3] 
except SiC and saturated drift velocity values. 
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Table 9 Calculated Figures of Merit at 300 K 

















APPENDIX B 


History and Nomenclature of Sic 
B . l Introduction 

Silicon carbide is an unusual semiconductor which 
rarely occurs in nature. In the solid state, it is the only 
known compound of silicon and carbon. Its existence has 
been known since the 19th century; it has yielded its 
secrets very slowly, however. 

This appendix will give a brief history of the 
metamorphosis of silicon carbide from something that 
occurred by chance in single crystal form as a byproduct of 
the Acheson process of making abrasive, to something that is 
now produced, doped or undoped, in the form of sublimation 
grown boules, using the modified Lely method, in continually 
increasing diameters. Currently, there are enough 
applications of Sic, so that once the wafer diameter exceeds 
three inches of usable surface area the material should be 
a competitive semiconductor. 

This appendix will also attempt to explain the 
nomenclature associated with silicon carbide, i.e. such 
things as poly type naming. There are in excess of 130 
different polytypes of silicon carbide that have been found 
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thus far; this is a small fraction of those that are 
theoretically possible. 


B.2 History 


The hardness of silicon carbide is between diamond and 
topaz [41] or sapphire [40] and only diamond and boron 
nitride are harder. 

The commercial value of silicon carbide [41] as a 
cutting and polishing abrasive was realized by A.H. and E.H. 
Cowless who had a U.S. patent for the production of silicon 
carbide in 1885 and A.G. Acheson who had a British patent in 
1892. 



UNREACTED 
MIXTURE 


CARBON 

CORE 



Figure 28 Acheson Furnace 


The Acheson process involves filling a rectangular 
volume, as seen in Figure 28, with the reactants: sand 50%, 
coke 40%, sawdust 7%, and NaCl 3%. Through the center of 
this mixture runs a carbon core which is the electrical 
heating element. The time -temperature heating cycle exceeds 
30 hours with a maximum temperature of 2700 °C. 
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During this process the sawdust shrivels, keeping the 
mixture porous. The salt reacts with impurities and the 
volatile reactants escape with the large volume of carbon 
monoxide generated, helping to improve the purity of the 
mixture. 

The mixture settles during the heating process forming 
voids in the Sic mass. The Acheson crystals are found in 
these voids. These crystals usually have an area of a few 
cm 2 with a thickness of a few mm. They usually have one 
well developed face. These crystals are heavily doped with 
a mixture of dopants and polytypes. The reaction results in 
predominantly one polytype, 6H. There are a large number of 
imperfections, e.g. twinning, and screw dislocations, in the 
crystal structure as well. 

These crystals were all that were available to the 
scientific community until 1955 when J.A. Lely published his 
method of growing purer crystals of known doping. 

The Lely method emulates the sublimation and 
recrystallization that occurs in the voids formed during the 
Acheson process on a laboratory scale. In later runs of the 
Lely method, a porous graphite cylinder is surrounded with 
the purest, light green, technical grade polycrystalline 
silicon carbide available (Fig. 29) . The charge of 
polycrystalline Sic in the Lely chamber is heated to 2500°C. 
in an argon atmosphere to which dopants can be added, i.e. 
dopants can be introduced by adding them to either the 



72 


POLYCRYSTALLINE SiC POROUS GRAPHITE 



SINGLE CRYSTAL SiC GRAPHITE 


Figure 29 Lely Chamber 

polycrystalline silicon carbide or to the argon. The SiC 
diffuses through the graphite cylinder and recrystallizes as 
single crystals on the interior of the porous cylinder. 

In practice, the yield is very low and the nucleation 
is uncontrolled. This was a tremendous improvement over the 
Acheson process in spite of its drawbacks. The yield was 
predominantly of the 6H polytype followed by the 15R and the 
4H polytypes [15,40]. 

The publication of two papers [1,10] using a seed 
crystal solved the nucleation problem of the Lely method. 
The resulting single crystal grows both in thickness and in 
width from the seed crystal. The modified Lely method 
(Figure 30) uses a porous graphite cylinder with the seed at 
the cooler part of the temperature gradient 2200 to 2300 [7] 
or 2400°C [l] . 
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The starting material is the purest Sic available or it 
could be compounds that contain either silicon or carbon or 
both. Doping is accomplished through the dopants in the 
starting material or by adding them to the atmosphere in the 
growth chamber. The three essential parameters are the 
temperature which determines the polytype, the temperature 
gradient and the pressure which determine the transport 
velocity, i.e. the growth velocity. High pressure is used 
to retard the growth until the proper temperature has been 
reached. 

The grown boules are cut into slices with a diamond 
watering saw and the exterior, of polycrystalline material 
that forms during cool down, is removed by coring the center 
of the slice. The wafers are then polished. Sic can be 
visualized as stacked layers, and the layers are double 
layers with a close packed layer of silicon and a layer of 
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carbon. Each of the carbon atoms is directly above each of 
the silicon atoms. Therefore, each wafer always has a 
carbon face and a silicon face. 

B.3 Structure and Nomenclature 

Silicon and carbon are both group four elements making 
SiC a IV-IV compound semiconductor. Each silicon atom is 
tetrahedrally surrounded by four carbon atoms and each 
carbon atom is likewise surrounded tetrahedrally by four 
silicon atoms. The atoms form parallel planes visualized as 
mentioned in the previous section, as double layers of 
silicon and carbon atoms when viewed along the 1120 plane 
[40,41]. 

If a plane of close packed spheres is viewed from above 
(or below) there are six spaces around each sphere and two 
ways to stack another plane of close packed spheres on the 
first plane; the second plane will occupy three of the six 
spaces. Using an A, B, C position notation, the first plane 
is the A position, the second plane is in the B position and 
the third can be in either the A or the C position (Fig. 
31) . Each additional plane can occupy one of the two 
positions not already occupied by the underlying plane, 
possibly forming large unit cells. The layers can form 
either cubic (trigonal) or hexagonal symmetry. As an 
example in cubic symmetry the adjacent layers of B are in 
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adjacent layers of B are in the same position (ABA) . 

The polytypes form either cubic (C) , also called beta- 
SiC, hexagonal (H) , or rhombohedral (R) structures; the 
latter two are also referred to as alpha-SiC. 

There are several different notations used to identify 
the crystal structure. The Ramsdell notation uses a number 
followed by a letter. The number represents the number of 
double layers in the unit cell and the letter represents the 
structure, e.g. cubic (C) , hexagonal (H) , or rhombohedral 
(R) . The Jagodz inski notation characterizes the layers by 
adjacent layers in the unit cell. That is whether they are 
cubic or hexagonal. 

Table 10 gives some examples of different polytypes 
(3C,6H) in different notations and the number of 
inequivalent sites. 
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Table 10 Poly type Notation 


RAMSDELL 

STACKING 

JAGODZ INSKI 

INEQ 

H 

C 





E 

U 





X 

B 

2H 

AB 

h(h) 

1 

1 

0 

3C 

ABC 

c(c) 

1 

0 


4H 

ABAC 

hc(c) 

2 


g 

6H 

ABCACB 

hcc(hcc) 

3 

D 

2 

15R 

ABCACACBCABACABCB 

hcchc(hcchc) 

5 

2 

3 


The bandgap extremes are 2.2 eV for 3C and 3.2 eV for 
2H. The bandgap is approximately inversely proportional to 
the fraction of each inequivalent site type. As an example, 
the 6H polytype has 1/3 of its inequivalent sites of the 
hexagonal type and 2/3 of its inequivalent sites of the 
cubic type, so 1/3 of 3.2 + 2/3 of 2.2 = 2.867. The bandgap 
of 6H is reported to be 2.93 eV. 































APPENDIX C 


Least Squares 


C.l introduction 

Designing and conducting an experiment to obtain the 
best possible data requires a significant investment of time 
and effort. The data then needs to be evaluated to extract 
results that are unique, unambiguous, and with the smallest 
possible error. The least squares approach yields values 
that are normally distributed about the true values with the 
smallest variance. This appendix will briefly explain the 
least squares method of analysis and how it was applied in 
this study. 

C. 2 Least Squares Method 

Least squares is a method of analysis, developed by 
Carl Friedrick Gauss in the 18th century (c.1795), where a 
set of parameters (Pj) are estimated from measurements of 
some other related quantities (Y { ) so that the difference 
between the experimental (or observed) values and the 
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adjusted (or calculated) values are minimized. 
Theoretically the adjusted values are closer to the true 
values than the observed values. 

The basis for the least squares method is that the 
errors between the observed data values and the calculated 
values can be described by the normal (Gaussian) error 
distribution function. When this is true the likelihood 
function takes the following form: 

l(P) - constant - S / (2a 2 ) (54) 

where 1(P) is the likelihood that the parameter P is the 
correct value, a 2 is the variance of the data and S is the 
sum of squares function given by 

S = (Yj - Y{X i ,P 1 ,P 2 Pj)) 2 (55) 


In the sum of squares function Y. is the observed value and 
the function Y (X { , P 1 , P 2 . . . Pj) is the calculated or adjusted 
value based on an appropriate modeling equation. 

Then the maximum likelihood function (dl(P) /dP] reduces 
to the method of least squares. This is the reason the 
weighting factor, the reciprocal of the variance, needs to 
be included in every calculation. 

The weighted sum of squares S is the fundamental 
equation of the least squares method. 



79 


5 = K-i (-4 [Yi ~ Y(X if P lt P 2 . . .Pj)] 2 ) (56) 

The object is to make S as small as possible. When S is 
minimized, by definition, the derivatives of S are all equal 
to zero. Then by taking the partial derivatives of S with 
respect to each parameter (dS/dPj) , a set of equations is 
produced which can be used to solve for the parameters of 
interest. 

C.3 Least Squares Applied 

This study follows the procedures outlined by Philip R. 
Bevington for the linearization of the fitting function in 
reference [53] to evaluate the data. He uses the gradient- 
expansion algorithm of Marquardt to combine the best 
features of the gradient search with the method of 
linearizing the fitting function. This is necessary since 
the gradient search is best when approaching the minimum 
from a distance, but converges slowly when near the minimum; 
and the analytical method of using a Taylor's expansion of 
the fitting function y(x) can only be used reliably when in 
the immediate vicinity of the minimum. The algorithm offers 
the advantage that only the simpler first-order expansion 
needs to be valid when near the minimum. 
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The algorithm increases the diagonal terms of the 
curvature matrix a which he defines as 


~ r l dYjXj dYUJ 
jk " a\ dPj dp k 


(57) 


by a factor A which determines which method predominates 
during an iteration. If X is very small, the solution is 
similar to a Taylor's expansion like the analytical method. 
If X is very large, the diagonal terms dominate giving n 
separate equations and the solution is similar to the 
gradient search method. 

Marquardt's recipe as given in ref. [53] is the 
following: 

1. Compute S(P). 

2. Start initially with X = 0.001. 

3. Compute £P and S(P+5P) with this choice of A,. 

4. If S (P+5P) > S (P) , increase A by a factor of 
10 and repeat step (3) . 

5. If S (P+5P) < S (P) , decrease A by a factor of 
10, consider P' = P+5P to be the new starting 
point, and return to step (3) substituting P' 
for P. 

This is done so that S decreases and A is small enough 
to make use of the analytical method. This algorithm was 
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implemented by the computer subroutine CURFIT which was 
modified to make it run more efficiently and without 
overflows 16 . 


16 The modified subroutine CURFIT was made available 
to me courtesy of Dr. Edward Haugland of NASA LeRC. 

C -<£ . 
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C SUBROUTINE CURFIT — FROM BEVINGTON'S BOOK, PG. 237 
C 

C PURPOSE 

C MAKE A LEAST-SQUARES F IT TO A NON-LINEAR FUNCTION 
C WITH A LINEARIZATION OF THE FITTING FUNCTION 

C 

C USAGE 

C CALL CURFIT (X, Y, SIGMAY, NPTS , NTERMS , MODE, A, 

C DELTAA , S IGMAA , FLAMDA, YF IT , CHISQR) 

C 

C DESCRIPTION OF PARAMETERS 

C X - ARRAY OF DATA POINTS FOR INDEPENDENT VARIABLE 

C Y - ARRAY OF DATA POINTS FOR DEPENDENT VARIABLE 

C SIGMAY- ARRAY OF STANDARD DEVIATIONS FOR Y DATA POINTS 
C NPTS - NUMBER OF PAIRS OF DATA POINTS 
C NTERMS- NUMBER OF PARAMETERS 

C MODE - DETERMINES THE METHOD OF WEIGHTING THE FIT 
C +1 (INSTRUMENTAL) WEIGHT (I) = 1/S IGMA ( I ) * * 2 

C 0 (NO WEIGHTING) WEIGHT (I) = 1 

C -1 (STATISTICAL) WEIGHT(I) = 1/Y(I) 

C A - ARRAY OF PARAMETERS 

C DELTAA- ARRAY OF INCREMENTS FOR PARAMETERS A 
C SIGMAA- ARRAY OF STANDARD DEVIATIONS FOR PARAMETERS A 
C FLAMDA- PROPORTION OF GRADIENT SEARCH INCLUDED 
C YF IT - ARRAY OF CALCULATED VALUES OF Y 
C CHISQR- REDUCED CHI SQUARE FOR FIT 
C 

C SUBROUTINES AND FUNCTION SUBPROGRAMS REQUIRED 
C FUNCTN (X , I , A ) 

C EVALUATES THE FITTING FUNCTION FOR THE ITH TERM 

C FCHISQ (Y, SIGMAY, NPTS, NFREE, MODE, YFIT) 

C EVALUATES REDUCED CHI SQUARE FOR FIT TO DATA 

C F DERI V (X, I , A, DELTAA, NTERMS, DERIV) 

C EVALUATES THE DERIVATIVES OF THE FITTING FUNCTION 

C FOR THE ITH TERM WITH RESPECT TO EACH PARAMETER 

C MATINV (ARRAY, NTERMS, DET) 

C INVERTS A SYMMETRIC TWO-DIMENSIONAL MATRIX OF DEGREE 

C NTERMS AND CALCULATES ITS DETERMINANT 

C 

C COMMENTS 

C DIMENSION STATEMENT VALID FOR NTERMS UP TO 10 
C SET FLAMDA = 0.001 AT BEGINNING OF SEARCH 

C CALCULATED VALUES OF SIGMAA (J) DEPEND ON VALUE OF FLAMDA 
C SET FLAMDA = 0.0 AFTER CONVERGENCE TO CALCULATE SIGMAA (J) 
C 

C CHANGES FROM BEVINGTON ORIGINAL 

C CHANGED STMTS 73 AND 84 SLIGHTLY TO PREVENT OVERFLOWS 
C (10/87) 
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C ADDED POSSIB ILTIY OF ENTERING F LAMDA = 0. (2/88) 

C ADDED STMT 28 TO ALLOW CALC OF WEIGHT (I) WHEN 
C S IGMAY (I ) =0 . (3/90) 

C INCREASED DIMENSION OF WEIGHT FROM 100 TO 300 (3/90) 

C ADDED PROVISION FOR ALPHA (J , K ) =0 . IN DO 73 LOOP (3/90) 

C ADDED PROVISION FOR TYPING VALUES OF ALL MATRIX VALUES 

C 

C EDITED 3/20/90 

C 

SUBROUTINE CURF IT (X , Y , S I GMAY , NPTS , NTERMS ,MODE , A , DELTA A , 
1 S IGMAA ,F LAMDA , YF IT , CHI SQR ) 

DOUBLE PRECISION ARRAY 

DIMENSION X (1) , Y (1) , S IGMAY ( 1 ) ,A (1) , DELTAA ( 1 ) ,S IGMAA (1) , 

1 YF IT (1) 

DIMENSION WEIGHT ( 300 ) , ALPHA ( 10 f 10) f B ETA ( 10 ) ,DERIV(10) , 



1 

ARRAY (10, 10) ,B (10) 

11 


NF REE =NPTS -NTERMS 

IF (NFREE) 13, 13,20 

13 


CHISQR=0. 

GO 

TO 

110 

C 



C EVALUATE WEIGHTS 

C 



20 


DO 30 11=1, NPTS 

1 = 

II 


21 


IF (MODE) 22,27,29 

22 


IF (Y (I) ) 25,27,23 

23 


WEIGHT (I ) -l./Y (I ) 

GO 

TO 

30 

25 


WEIGHT ( I ) =1 . / (-Y (I) ) 

GO 

TO 

30 

27 


WEIGHT (I ) =1. 

GO 

TO 

30 

28 


IF (SIGMAY (I ) ) 29,27,29 

29 


WEIGHT (I ) =1. /SIGMAY (I ) **2 

30 


CONTINUE 


C 

C EVALUATE ALPHA AND BETA MATRICES 
C 

31 DO 34 J = l, NTERMS 

BETA ( J) =0. 

DO 34 K=1 , J 

34 ALPHA (J,K) =0. 

41 DO 50 11=1, NPTS 

I = II 

CALL F DERI V (X , I , A , DELTAA , NTERMS , DERIV) 

DO 46 J = 1 , NTERMS 

BETA (J) =B ETA ( J) +WEIGHT (I ) * (Y (I ) -FUNCTN (X, I , A) ) *DERIV (J) 



DO 46 K=1 , J 

ALPHA ( J , K ) =ALPHA (J,K)+WEIGHT (I ) *DERIV (J)*DERIV (K) 

46 CONTINUE 

50 CONTINUE 

51 DO 53 J = 1 , NTERMS 
DO 53 K=1 , J 

53 ALPHA (K,J)=ALPHA(J,K) 

C 

C TYPE OUT ALPHA MATRIX (IF DESIRED) 

C (FORMAT FOR NTERMS UP TO 7) 

C 

C TYPE 550 
C DO 56 J = l, NTERMS 

C 56 TYPE 560, (ALPHA (J,K) ,K = 1, NTERMS) 

C TYPE 550 
C 550 FORMAT (/) 

C 560 FORMAT (5X, 1PE11. 3, 6E11. 3) 

C 

C EVALUATE CHI SQUARE AT STARTING POINT 
C 

61 DO 62 11=1, NPTS 
I=II 

62 YF IT (I ) =FUNCTN (X, I ,A) 

63 CHISQl=FCHISQ ( Y , S IGMAY , NPTS , NF REE , MODE , YF IT) 
C 

C EVALUATE MODIFIED CURVATURE MATRIX ARRAY 
C 

71 DO 74 J=l, NTERMS 
DO 73 K=l, NTERMS 

AX=SQRT (ALPHA (J, J) ) *SQRT (ALPHA (K, K) ) 

IF (AX) 70,70,72 
70 ARRAY (J,K) =0. 

LOOP=73 

TYPE 200, J,J,K,K, LOOP 
GOTO 73 

72 ARRAY ( J , K) = ALPHA ( J , K) /AX 

73 CONTINUE 

74 ARRAY ( J , J ) =1 . +F LAMDA 
C 

C TYPE MODIFIED CURVATURE MATRIX, IF DESIRED 
C 

C TYPE 550 
C DO 76 J=l, NTERMS 

C 76 TYPE 560, (ARRAY (J,K) ,K=1, NTERMS) 

C 

C INVERT MODIFIED CURVATURE (ERROR) 

C MATRIX TO FIND NEW PARAMETERS 
C 
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80 CALL MATINV (ARRAY, NTERMS, DET) 

C 

C TYPE MODIFIED CURVATURE MATRIX, IF DESIRED 
C 

C TYPE 550 
C DO 86 J=1 , NTERMS 

C 86 TYPE 560, (ARRAY (J,K) ,K = 1, NTERMS) 

C 

81 DO 84 J=l, NTERMS 

B ( J ) =A ( J ) 

DO 84 K=l, NTERMS 

AX=SQRT (ALPHA (J, J) ) *SQRT (ALPHA (K, K) ) 

IF (AX) 82,82,83 

82 LOOP=84 

TYPE 200, J,J,K,K, LOOP 
GOTO 84 

8 3 B (J) =B ( J ) +B ETA (K ) *ARRA Y (J,K)/AX 

84 CONTINUE 

C 

C IF CHI SQUARE INCREASED, INCREASE F LAMDA AND TRY AGAIN 
C 

91 DO 92 11=1, NPTS 
I =1 1 

92 YF IT ( I ) =F UNCTN (X , I ,B ) 

93 CHISQR=FCHISQ ( Y , S IGMAY , NPTS , NF REE , MODE , YF IT ) 

IF (F LAMDA) 94,101,94 

94 IF (CHISQ1-CHISQR) 95,101,101 

95 F LAMDA =10. *F LAMDA 
GO TO 71 

C 

C EVALUATE PARAMETERS AND UNCERTAINTIES USING ERROR MATRIX 
C 

101 DO 103 J=l, NTERMS 

LOOP=103 
A ( J ) =B ( J ) 

IF (ARRAY (J,J) .LE. 0. ) TYPE 201 , J , J , ARRAY (J , J ), LOOP 
IF (ALPHA (J,J) .LE. 0. ) TYPE 202 , J , J , ALPHA (J , J ), LOOP 
IF (ALPHA (J, J) .EQ. 0. ) ALPHA (J, J)=. 1 
103 SIGMAA ( J ) =SQRT (ABS (ARRAY (J, J) )/ABS (ALPHA (J, J) ) ) 

IF (F LAMDA) 104,110,104 


104 

F LAMDA =F LAMDA/10 . 


110 

RETURN 


200 

FORMAT (T2 'ALPHA ('11,11') 

OR ALPHA ('11,11') IS 

ZERO! DO 

LOOP' ,14) 

201 

FORMAT (T 2 'ARRAY ('11,11') 

= ' , 1PE10 . 3 , ' DO LOOP', 14) 

202 

END 

FORMAT (T2 'ALPHA ('11,11') 

= ' , 1PE 10 . 3 , ' DO LOOP', 14) 



APPENDIX D 


Cryostat Schematics 



Figure 32 High temperature furnace and sample 
holder. 
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TOP VIEW NOT TO SCALE 
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Figure 33 Firebrick machining dimensions 
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NOT TO SCALE 
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Figure 36 Top of the sample platform and thermocouple 
support . 





APPENDIX E 


Effective Mass 


The density of states function is used to calculate the 
carrier concentrations in the material. To do this the 
density of states effective mass needs to be known. 

There are two ways of reporting the effective mass in 
the literature. One group [6,18,20,21,34] reports results 
that are consistent with 


Ml . ( 


m , 


2 

m t 


)l/3 


( 58 ) 


where m t is the longitudinal effective mass, m t is the 
transverse effective mass, m* is the effective mass of 
electrons or holes per conduction band minimum or valence 
band maximum, respectively. The other group [7,30,32,33,35] 
reports results that are consistent with [36], 


ir>d 


= M 2/ 3 [;i7j*jnt] 1/3 


( 59 ) 


where M is the number of conduction band minima or valence 
band maximum for electrons or holes, respectively, and m d * 
gives the total density of states effective mass for 
electrons or holes. 
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These are used in the calculation of the effective 
density of states function N c v given by 

N c> v = 2M [ iMlJSl ] 3/z ( 60 ) 


Here M is included if equation (58) is used and dropped if 
equation (59) is used for the effective mass, k is 
Boltzmann's constant, T is the temperature in kelvins, and 
h is Planck's constant. The subscripts c or v stand for the 
conduction or valence band respectively. 

H. J. van Daal [8] also made a plot, effectively, of 


( M * mg 
9 


( 61 ) 


versus the reciprocal distance between minority centers. 
Here, M is the number of valleys in the valance band, m^* is 
the density of states effective mass of holes, and g is the 
degeneracy factor. For aluminum, he found (M*m d */g) =0 . 5 
which, for M=l, m d *=l, would give g=2. He, also, found the 
effective mass of electrons to range from 0.72m to 1.0m . 

In conjunction with this, there is a controversy in the 
SiC community concerning the number of conduction band 
minima that are within the first Brillouin zone. There 
appear to be six conduction band minima in Sic, just as 
there are in silicon, but are there 3 or 6 valleys within 
the first Brillouin zone [20,30,35]? As an example, 
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germanium [36] has eight conduction band minima but only 
four are within the first Brillouin zone. Therefore, four 
is used for M in the effective mass or density of states 
equation for germanium. 

The uncertainties in the number of conduction band 
minima being 3 or 6, and in the magnitudes of the transverse 
(0.24 to 0 . 3 5 ) m 0 and longitudinal (0.34 to 1.5)m o effective 
masses [21,22,30,35,37,38], could be unsettling. But 
fortunately, when the total density of states effective mass 
is calculated using these diverse quantities there is 
agreement that m^* = (1.0±0.2)m o for electrons in the 
conduction band; either directly by the reported normalized 
effective mass [7,9,30,32,33,35] or by calculating m de * using 
equation (59) and a three valley approach using the reported 
per valley effective mass [6,18,20,21,34]. In this research 
m* = 0.44 m was used after Muench et al [6]. Therefore m,* 
= M z/3 (0.44 m 0 ) = 0.915 m 0 with M=3. Since the valence band 
maxima is at the center of the Brillouin zone [15,22,23] m dh 
= (1.0)m o [8,9,23,34,39]. 
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